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Free Electron Laser 



This paper treats the growth of the radiation field in two-dimensions of a free electron 
laser on an electron beam with finite axial and transverse dimensions in the amplifying configuration. 
The general, self-consistent, nonlinear analysis includes various efficiency enhancement schemes, 
diffraction and refraction. In the axially symmetric, low gain, resonant macro particle limit, we ob¬ 
tain an analytical expression for the gain. An illustration at 10.6 ^ is given. 
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2-D Nonlinear Theory of the Free 
Electron Laser Amplifier for an 
Electron Beam with Finite 
Axial and Transverse Dimensions 

Many current experiments of the free electron laser (FED, ultilize electron beams from a 
millimeter to a few centimeters in pulse length. The short pulse length is typical of high energy 
accelerators such as RF Linacs and microtrons. The finite length effect of the electron beam on 
the radiation was found to be important in the Stanford oscillator experiment.^” Currently, 
many experiments in the amplifying mode are being conducted with the short electron beam 
pulses. We have an analytical expression for the gain pulse of the radiation field, applicable to 
these experiments. 

One-dimensional analysis of the radiation field for electron beams of finite length have 
been numerically simulated on computers.^*"*’ The effect of the finite transverse dimensions 
was either not included, or incorporated through filling factors. On the other hand, previous 
three-dimensions self-consistent formulation^” of the radiation field for a semi-infinitely long 
electron beam in the amplifying configuration has resulted in a number of interesting effects 
not obtainable by the 1-D formulation. Numerical effort to find the growth of the 3-D radia¬ 
tion field on the finite length electron beam began only recently. 

In this paper, we will present a fully 2-D, self-consistent, non-linear, analytical analysis of 
the FEL process in the amplifier mode of operation treating the finite length and transverse 
effects associated with both the electron beam and the radiation beam. Our formulism also 
includes various efficiency enhancement schemes: (i) contouring in the longitudinal direction 
the amplitude and/or the wavelength of the magnetic wiggler field, and (ii) applying an external 
MaauKfipt nbinitted January 25, 1982. 
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D.C. electric field. Analytical results in the amplifying configuration are obtained in the low 
gain, trapped particle regime. 

The schematic of the configuration is shown in Fig. I. The generalized vector potentials 
of the right-handed, helicc!, static magnetic wiggler field and the electromagnetic radiation field 
are 

A» (z) - (z) [cos (J*^ (z') dz') e„ + sin ( k^. (z') dz’) ej (\) 

Ar (jf. y. z. /) - An (x. y, z. t) [cos (— z - at t A <p ix. y. z. i)) 

- sin (— z - <0 r + ^ (x y. z. /)) e,I (2) 

c 

where A^ and are all slowly varying amplitude and wave number of the wiggler field and Ag 
and If are slowly varying amplitude and phase of the electromagnetic radiation field following 
the electron pulse. We also include an external DC electric field, K 

for the purpose of efficiency enhancement. 

In this analysis we will not consider the gradient in the wiggler Held. This is a good 
approximation if k^r^ « f, where r* is the radius of the electron beam. If the FEL is operat¬ 
ing in a trapped particle mode, we also require**’ < (-y^o (8 V2 

where y,o ” (1 - - If I/4»/(yomo<^^). yo “ yzoyq- ya “ 

(1 + If 1^ /4,^(0)/(mo and v,© is the axial velocity at z * 0. 

The electron motion can be described in terms of their phase ifr in the ponderomotive 

wave: 



where ^(xo, yo- fo- ” Jl ®' {k^(z’) A- «/c)<fe' -F <ur -I- -Pfx, y, , z, t) is the phase for 

the electron, which was at (xq, y^, ^o) at / — 0, y-y,yi, y, •• (1 - v//fy^ - 
(1 + If P/4»(z)/(m^c^))'^^ V, - [d^ldt- + o)l/U*(z) + w/cl is the axial velocity, z — 
2 




NRL MEMORANDUM REPORT 4774 

^0 + /j V, (jfQ, >' 0 . ^ 0 . is the axial position of the electron and fo is the axial position of 
the electron relative to the center of the electron beam at / — 0 . 


The wave equation for the radiation field is (V^ - c ^dVdt^)Ag - -4flrc"'J , where 
/— yo)h((o^ 

-y“‘Pi8(x - Xo)S(>'- >o)S(z - z) (4) 

is the current, (xq. >> 0 ) are the particle's transverse positions at r - 0 , 0 (xq. >> 0 ) is the 
transverse current density profile, /i(lo) 's the macroscopic electron pulse shape, no is the peak 

iff I 

current density, and Pi = is the transverse momentum. 


We can rewrite the radiation field as A« — a^fx, y, z, t) exp [/(uz/c - o»t)le+ + c.c., 
where a, •- Aji exp(/<p) is the complex amplitude of the radiation field, and e± - (e^ ± &,)I2 
is a new coordinate system. The wave equation assumes the form 


dyt c |dz c dr 




~ f“ 


exp j-i| (A:„(z')+ <i»/c)<fe'- w^jsfx- Xo)80'- ^' 0 ) 6 ( 2 - f). (5) 

In obtaining (5), we have used the fact that KdVdz^ - c"^dV9f^)a«l « WcKd/dz + 
c-'d/ar)fl*l. 


The solution for an can be separated into the input radiation field a„, and the excited 
radiation field, a„, such that an — a,„ + a„- The excited radiation field can be shown to be 

- L /_«‘too >- 0 . fo. X. y. z, r')6(z - c(f- r') - z) (6) / 

where / 

(k»(r') + «/c) dr' - ^ 


f “ 0 (jfo- >o)d (fo) ■ 


(r - f')y 

IJ (x-xq)^ + 0'-yo)^ | « 

‘*11 2 c(r-r’) j c 


Hr 
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The integral in time of Eq. (6) can be evaluated by changing the argument of the delta func¬ 
tion. 


■ /o * /-=. ‘*^0 /-»J'o. fo. X. y. z. /') — - — 

di' 




where to(xo. yo, (q, i) is the retarded time associated with the electron, which originated at 
(jfo- > 0 . ^ 0 ^ ^ ” 0- The retarded time satisfies the equation 


^0 + Zr^’’o) + CU - Tq) - ^ -H Zc(t) (8) 

where \„U')dt' is the macroscopic location of the center of the electron beam at 

time I, v^(t) «• ot/(k„(Zc(t)) + <u/c) is the macroscopic velocity of the electron pulse, f is the 
position of the electron relative to the center of the electron beam at time t. 


The complex radiation amplitude in Eq. (8) can be evaluated if we make the following 
simplifying assumptions. For experimental parameters of interest, we can assume that the 
bunching mechanism does not alter the macroscopic electron pulse shape, hence, it travels 
undistorted through the interaction region. We will assume that the electron beam has an axi¬ 
ally symmetric Gaussian profile in the transverse direction, i.e., 9 (xq, yo) “ 
exp[-(xo^ +yo )/ro^). Furthermore, we will assume that the waist of the input radiation field tq 
is much larger than the radius of the electron beam r*, such that i is approximately a function 
of $0 and t only. The excited radiation field takes the form 

^ AJi + z,U)- cO-Tq)) (l-FVc)y/ 

y c(r-To)-K* 

«p[' To) + ^ (r. f. ro))l 
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where zt, - ri m/lc is the Rayleigh length associated with the electron beam radius. Equations 
(3) and (9) describe self-consistently a general, nonlinear, 2-D, PEL amplifier with a macros¬ 
copic pulse shape 

For the purpose of illustrating the finite length pulse efifects in an PEL amplifier operating 
in the low gain limit, i.e., |ni„| » loexi^ 've take the electron beam profile to be uniform, i.e., 
hdo) “ 1 for l^ol ^ f-»/2 and “ 0 for ||ol > f-*/2, where is the length of the elec¬ 
tron pulse. We also make the constant phase, resonant particle approximation. In this approxi¬ 
mation all particles are assumed to have the same constant phase, The electron beam in 
this approximation consists of a pulse train of macro particles separated in distance by 2n\Qj<,t. 
Furthermore, we will limit ourselves at this point to a constant parameter wiggler and consider 
only an external DC electric potential. The amplitude and phase of the total field are 

t) - Ai„ - cosin + 4 sin 0;^) (10a) 

cos^K ~ h sin(10b) 
where 4,^ — /, - Ke (I), I, - Im (/). / - £jj““q/| - Ei is the 

exponential integral function, 

qi - (-(z»)(2>/fi(f - fo.,) - K*)"'. 

9i/ ” (-»*)(2y//?(f - fo.i,) -c*)"'. 

fo.K ■“ ( (for f < f'*/2) and i - L*/2 (forf > LJ2) is the upper limit of the integration, f o./ 

- f -(c- v„)r (for ( -(c- v„)/ > -L^/2) and (o.t - - V2 (for f -(c- v„)/ < -Lt/2) is 
the lower limit of the integration, and is the resonant gamma associated with the axial 
motion. 

A more realistic electron beam profile h({o) "* f -(2{(/Lfc)^ (for fo ^ EiJ2) and h({o) 

- 0 (for Ifol ^ f'»/2) can also be integrated. The result is not given here, because the more 
complicated expressions would obstruct the initial understanding of the physical process of the 
pulse propagation. 
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As an example of a 10.6 /im FEL utilizing a CO 2 laser as an input field, we choose an 
electron beam of energy 25 MeV (yo — 50), current of / — 5 A and radius (Gaussian profile) 
of th — 0.5 mm and pulse length = 3 mm. Such a beam has a peak density on axis of /iq 
- 1.3 X 10 “ cm“^ (ait — 2.0 X 10’®sec“'). The constant parameter wiggler has a magnitude 
of - 5.0 kG and wavelength of /* = 2.8 cm which gives - 2.2 x 10^ statvolts. The 
wiggler velocity is voi " 2.6 x 10“^ c which gives yi “ 1.35 and y^ — 37. The input CO 2 
power density is taken to the P,, -=4x10* W/cm^ which gives - 0.30 statvolts. Our illus¬ 
tration assumes resonant macro particle approximation and an applied D.C. electric potential 
such that sin 0 ft = 0.6. 

The schematics of the gain 

G(r.(.t)~ (AK(r.(. t) - Ai,)/A,„ 

are shown in Figs. 2 and 3. The slashed bars in the (2, t) plot of Fig. 2 denote the locations of 
the electron beams at tj - 1 m/c and r — 2 m/c, which c is the speed of light. The solid lines 
in the (r, 1) plot are the light lines. The gain pulse on axis are plotted at times r, and f 2 - We 
see that the excited radiation pulse grows and spreads beyond the electron beam pulse. The 
transverse variation of the gain at f — 0 for various times are plotted in Fig. 3. The decrease 
of radiation field far from the axis is due to refraction toward the center of the beam. 

We have obtained a general expression for the growth of the 2-D, stimulated radiation 
pulse on an electron beam of finite axial and transverse dimensions in an FEL amplifier. We 
included diffraction as well as refraction. In the axially symmetric, low gain, resonant macro 
particle limit, we have an analytical expression for the radiation gain. The formalism presented 
here can be modified to study the radiation build up and "laser lethargy” in the FEL oscillator. 
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